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Abstract 

In this work we consider a size-structured cannibalism model with the 
model ingredients (fertility, growth, and mortality rate) depending on size 
(ranging over an infinite domain) and on a general function of the stand- 
ing population (environmental feedback). Our focus is on the asymptotic 
behavior of the system, in particular on the effect of cannibalism on the 
long-term dynamics. To this end, we formally linearize the system about 
steady state and establish conditions in terms of the model ingredients 
which yield uniform exponential stability of the governing linear semi- 
group. We also show how the point spectrum of the linearized semigroup 
generator can be characterized in the special case of a separable attack 
rate and establish a general instability result. Further spectral analysis 
allows us to give conditions for asynchronous exponential growth of the 
linear semigroup. 

Keywords: Size-structured populations; cannibalism; linear semigroup 
methods; asynchronous exponential growth 



1. Introduction 

Cannibalism is a phenomenon observed among many species, e.g. certain fish 
populations. Sophisticated population models are capable to elucidate a poten- 
tially stabilizing effect of cannibalism, underscoring that certain populations 
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may benefit from cannibalism when resources are limited. Consequently, the 
effects of cannibalism on the long-term dynamics of populations have attracted 
considerable interest and have been analyzed for various structured population 
models (see [5j[T4] for further references). 

Structured population models are typically formulated as partial differential 
equations for population densities. Diekmann et al. have developed a general 
mathematical framework to study analytical questions for structured popula- 
tions (see [3 13 |9]), including those pertaining to linear /nonlinear stability of 
population equilibria. In this context it was recently proven for large classes of 
structured population models, formulated as integral (or delay) equations, that 
the nonlinear stability/instability of a population equilibrium is completely de- 
termined by its linear stability/instability, a result commonly referred to as the 
"Principle of Linearized Stability" . 

In a series of recent papers we have successfully applied linear semigroup 
and spectral methods to formulate biologically interpretable conditions for the 
linear stability/instability of equilibria of various structured population models, 
see [TTJ fTS] and the references therein. In these problems we assumed that 
any effect of intraspecific competition between individuals of different sizes on 
individual behavior is primarily due to a change in population size and that every 
individual in the population can influence the vital rates of other individuals 
("scramble competition"). 

When the competition among individuals is based upon some hierarchy in 
the population, often related to the size of individuals, environmental feedback 
is incorporated in the model through infinite dimensional interaction variables 
( "contest competition" ) . This situation in its simplest form is of relevance for a 
forest consisting of tree individuals in which the height of a tree determines its 
rank in the population. Taller individuals have higher efficiency when competing 
for resources such as light, while individuals of lower rank do not influence the 
vital rates of individuals of higher rank. 

Of interest in this work is the linear stability analysis of population equilibria 
of a continuous quasilinear size-structured model, recently discussed in [7]. In 
particular, it was proven there that the Principle of Linearized Stability holds 
true for this model. The density evolution of individuals of size s is assumed to 
be governed by the partial differential equation 



n t (s, t) + ( 7 (s, E(s, t j) n(s, t)) s + (/x(s, E{s, t)) + M(s, t)) n{s, t) = 0, (1.1) 



defined for s G (0, oo) and t > 0. The density of zero (or minimal) size individ- 
uals is given by the nonlocal boundary condition 




t > 0. 



(1.2) 



The initial condition reads 



n(s, 0) = uq(s), s e [0, oo). 



(1.3) 
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Here /3, \i and 7 denote the fertility, mortality and growth rate of individuals, 
respectively. All vital rates are size-dependent. Moreover, it is assumed that 
the mortality and growth rates of individuals depend on the function 



Here a denotes the size-specific attack rate, while c denotes the energetic value 
of the attacked prey. Hence the function E models the assumption that the extra 
energy available due to cannibalism is channeled into the growth of individuals 
and affects their starvation-driven mortality. The extra size-specific mortality 
rate due to cannibalism is given by 



(Note the switch of variables in a in contrast to (II. 4|) .) We remark that the can- 
nibalism model discussed here incorporates the environmental feedback variable 
used in the hierarchical size-structured model in [13] . 

We impose the following regularity conditions on the model ingredients: 

• li = n{s, E) e C b ([0, 00); C£[0, 00)), n > 

• 7 — 7( s ; E) S Cb ([Oj 00); Cl [0, 00)), 7 > 70 > for some constant 70 

• P = p{s) eC 6 ([0,oo)), /3>0 

• c = c(s) G C 6 ([0,oo)), c > 

• a = a(y, s) e C b ([0, 00); ^([0, 00))), a > 0. 

The subscript b indicates that functions (and derivatives in case of C 1 ) are 
bounded. For clarity in later developments we will write D^a for the derivative 
of a with respect to its second argument. The regularity assumptions above 
are tailored toward the linear analysis of this work. They might, however, not 
suffice to guarantee the existence and uniqueness of solutions of Eqs. (|1.1[) - (|1.5|) . 
even in the steady-state case. Well-posedeness of structured partial differential 
equation models with infinite dimensional environmental feedback variables is 
in general an open question. It has recently been shown in [T] that population 
models with infinite dimensional interaction variables may exhibit a more com- 
plicated dynamical behavior than the simple size-structured model of scramble 
competition. 

The study of hierarchical models in the literature (see e.g. [3] and the ref- 
erences therein) is largely based on a decoupling of the total population quan- 
tity from the governing equations and a transformation of the nonlocal partial 
differential equation (jl . 1|) into a local one. This technique allows to prove well- 
posedness and to study the asymptotic behavior of solutions by means of ODE 
methods. For Eqs. (|1.1|) - (|1.5|) this approach fails, however, since the mortality 




(1.4) 




(1.5) 
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and growth rates depend on both size s and on the environment E(s, t). There- 
fore it seems unavoidable to study the original partial differential equation (jTTTJ) 
with the nonlocal integral boundary condition (jl.2D directly. Moreover, studying 
the linear stability of stationary solutions for models with infinite dimensional 
interaction variables by spectral analysis (as done previously for simpler cases 
in [TU 021 HO]) ^ as P roven difficult since eigenvalues are not given by an explic- 
itly available characteristic equation (see [13] ). Therefore we devised a different 
approach in [13] to establish a linear stability condition for a particular model 
with infinite dimensional interaction variables. 

Clearly, system (|1.1[) - (|1.5|) admits the trivial solution n = 0. In general, 
system (jl.l[) - (|1.5l) yields for a stationary solution n* 

n .(0) 7 (0, £.(0)) f f ii(y, E,(y)) + M,(y) \ 

where 



7(s,^(a)) r I J l{y,E*{y)) 



p oo poo 

E*(s)= I c(y)a(y,s)n s> ,(y)dy, M»(s) = / a(s,y)n*(y)dy. (1.7) 
Jo Jo 

Here and later on, starred quantities are stationary counterparts of the time- 
dependent functions in Eqs. (|l.ip - (|1.5p . For obvious reasons, we shall exclu- 
sively consider positive stationary solutions of the form (|1.6[) or the trivial 
solution n* = in the following. Moreover, to be consistent with later de- 
velopments, we shall assume throughout that stationary solutions n* have the 
regularity H /1:1 (0, oo). 

For fixed model ingredients (3, fj,, 7, a, c, a positive stationary solution n* 
with E* and M* as before satisfies the equation 

m r r^ E :^t^ dy ) ds =i. (1.8) 



o 7(s,£'*(s)) [ 7 j{y,E*(y)) 

Therefore, for solutions n = n(s,t) of the governing equations, it is natural to 
define the functional 

R{n)= - ( T-i / — TT ex P1-/ 7 — ^7 — u * (1-9) 

Jo 7{s,E{8,-)) { Jo j{y,E{y,-)) J 

where n, of course, determines the quantities E and M via (|1.4p . (1 1 . 5|) . i? may 
be regarded as the net reproduction rate of the standing population. Note that 
Eq. (|1.8[) requires ii(n*) = 1 for any positive stationary solution n*. 

In [9] a general framework was developed to establish the existence of steady 
states for general physiologically structured population models, however with 
finite dimensional interaction variables. For models with infinite dimensional 
interaction variables one can usually not formulate elegant necessary and suf- 
ficient conditions for the existence of steady state solutions. In this situation, 
one can construct positive stationary solutions by perturbation of solutions with 
a = const. In the latter case, the interaction variables E and M are constant, 
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hence the results in [S] apply. For an alternative approach we refer the reader 
to [2] . Throughout the rest of the paper we will tacitly assume that stationary 
solutions of the required regularity are available. 



2. Stability via dissipativity 

Given a stationary solution n* , we formally linearize the governing equations 
by introducing the infinitesimal perturbation u — u(s, t) and making the ansatz 
n = u + n*. After inserting this expression in the governing equations and 
omitting all nonlinear terms, we obtain the linearized problem 

u t (s,t) + (7* (s) u(s,t) + 7 _e(s, E*(s)) n*(s) F(s,t)) s + /i(s, E*(s)) u(s,t) 
+ (i B (s,E*(s))n*(s)F(s,t)+M*(s)u(s,t)+n*(s)N(s,t) = 0, (2.1) 



7*(0)«(0,t)= / ( J 9(a)-7K(0,S.(0))n.(0)c(«)o(a > 0))u(«,t)da, (2.2) 
Jo 

where we have set 

poo poo 

F(s,t)= c{y)a{y,s)u(y,t)dy, N(s,t) = a(s,y) u(y,t) dy, (2.3) 
Jo Jo 

and 

7*(s)=7(s, £*(*)), «G[0,oo). (2.4) 

We denote the Lebesgue space i 1 (0, 00) with its usual norm || ■ || by X and 
introduce the bounded linear functional A on X by 

Jo V7*(0) 7*(0) c(s)a( S ,0)J u(«)ds. (2.5) 

Next we define the operators 

Au = -7*(-)w«, Dom(^) = {tteW 1,:L (0,oo)|u(0) = A(u)}, (2.6) 
Bu =-0*0, £?.(.))+ 7^(0 + ^(0) «=-P%)« 011 ( 2 - 7 ) 



Cu=- / u(y) [c(y)a(y,.) (( 7b (., £*(.)) n*(.)) s + Mir(-, £*(•)) «♦(•)) 
Jo 

+c(y)D 2 a(y,.) 7E (;E*(.))n*{.)+a(.,y)n*(.)]dy onX. (2.8) 

Our regularity assumptions on the model ingredients and the stationary solution 
guarantee that these operators are well-defined, that the operators B and C are 
bounded on X, and that the operator A is closed and densely defined on X. 
Thus, the linearized system (j2.1jl . ()2.2() assumes the form of an initial value 
problem for an ordinary differential equation on X 

4«= (A + B + C) u, (2.9) 
at 
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with the initial condition 

u(0) = u . (2.10) 

Our objective in the following is to apply the Lumer-Phillips Theorem from 
linear semigroup theory (see 18J ) . To obtain stability by virtue of this result, we 
will extend our approach in [13 , which was previously devised for elongational 
flows in [15] . 

Theorem 2.1 The operator A + B + C is the infinitesimal generator of a quasi- 
contraction semigroup T = {T(t)} t >o of bounded linear operators on X. The 
semigroup is uniformly exponentially stable if 

n(s,E*(s))+M*{s) > \0{s) - 7JS (0,£,(0)) n„(0)c(s)a(s, 0)| 

\c{s) a(s, y) [(713(3/, E*(y)) n*(y)) y + fJ. E (y, E*(y)) n*(y)} 

+ c{s) D 2 a(s,y)-f E (y,E :t (y))n^(y) + a(y, s) n*(y)\ dy, s>0. (2.11) 

Proof. Assume that, for given h 6 X and fixed n S K, u G Dom(^4) is such 
that, for some A > 0, 

u-\{A+B + C + Kl)u = h. (2.12) 

Then we have 



u(s) sgnu(s) ds 

/•OO 

h(s) sgnu(s) ds — A / (j*(s) u(s)) s sgnu(s) ds 
Jo 

poo 

+ A / (k — n(s, E*(s)) — M*(s)) u(s) sgnu(s) ds 



Jo 

r-OC r-OO 

_A / / u {y)[c{y)oL{y 1 s){{'-iE{s,E sr (s))n^(s)) s + ^Eis.^E^s^n^s)) 
Jo Jo 

+ c{y) D 2 a(y, s) 7e(s, K*(s)) n*(s) + a(s, y) n*(s)] sgnw(s) ds. 

(2.13) 

First by changing the order of integration, we obtain 

pOO pOO 

-A / / u(y) [c(y)a(y, s) (('y E (s,E*(s))n*{s)) a + n E (s, E*(s)) n*{s)) 
Jo Jo 

+ c(y) D 2 a(y, s) Je(s, E*(s)) n*(s) + a(s, y) n*(s)} dy sgnu(s) ds 

p CO /*OC 

<A/ |u(s)| / [c{s)a{s,y) {{ lE {y-,E„{y))n tf (y)) y + H E (y, E*(y)) n*(y)) 
Jo Jo 

+ c(s) £> 2 «(s, y) 713(2/, S*(j/)) n*{y) + a(y, s) n*(y)] dyds. (2.14) 

Then we note that the set of points in the interval (0, 00) where u is nonzero 
is the countable union of disjoint open intervals (<2j,&t) with ai £ [0, 00) and 
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bi 6 (0, oo] such that on each of these intervals either u > or u < holds true 
with u{a,i) = unless a* = 0, and u(bi) = unless bi = oo. If (a*, bi) is any such 
interval on which u > we have 



h(s) sgnu(s) ds — A / (j*(s) u(s)) s sgnu(s) ds 



A / (k — fi(s, E*(s)) — M*(s)) u(s) sgnu(s)e?s 



/ \h(s)\ds + Xj (k- fj,(s,E*(s)) - M*(s))u(s)ds 
+ A 7* (a.,) it(aj). 
Similarly, on any interval (a*, 6^) where u < we have 

/i(s) sgn u(s) ds — A / (7 (s) u(s)) s sgnu(s) ds 

J at 

+ A / (/c — ^(s, E*(s)) — M*(s)) u(s) sgnu(s)ds 

< / |/i(s)| ds + A / (k - fi(s,E*(s)) - M*(s)) |u(s)|ds 
- A7*(ai) u(ai). 
Combining inequalities (|2.14p - (12.16[) we get the estimate 



||«||=5^ / u{s)sgnu(s)ds< ||/i|| + A7*(0)|tt(0)| 

/■ OO /"OG 

+ Ay \u(s)\ (k- n(s,E*(s))- M*(s)+ J \c(s)a(s,y) 

x [(7B(y,^*(y))"*(y))a + HE{y,E*{y)) n*(y)] 
+ c(s) D 2 a(s,y) j B (y, E*(y)) n*(y) + s) n*(y) | dy) ds 

Next we note that 

/3(s) 7E (0,£,(0K(0) 



(2.15) 



(2.16) 



(2.17) 



K0)| = |A( U )| < 



-c(s)a(s, 0) 



u(s)|ds. (2.18) 



7*(0) 7*(0) 
Now choose set such that, for s e [0, 00), 

K < fi(s,E*{s))+M*{s) - |/3(s) - lE (0,£,(0)) n»(0) c(a)a(s, 0)| 

/•OO 

- / |c(s)a(s,y) [(7^(2/, E*{y)) n*(y)) y + fi E (y, E*(y)) n*{y)] 
Jo 

+ c(s) D 2 a(s, y) j E (y, E*(y)) n*(y) + a(y, s) n»(y) I dy. (2.19) 
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For such k, we have the desired inequality 

[|u|| < for A > 0, (2.20) 

thus establishing dissipativity. We observe that the operator A + B + C + kI 
is densely defined and that the equation 

(XI-A)u = f (2.21) 

for f € X and A > sufficiently large has a unique solution u G Dom.4, given 

by 

lt(s)=expl — f — — dy 



with 



o i*{y) 

*^ )+ LAL $F) ir] (2 - 22) 



A( " )= ( 1 - A ( DP {-Xrro* 

A (Mjf^*-I?5y*}$H- (2 - 23) 

The fact that u € DomA is well defined by (|2.22[) . (|2.23p follows immediately 
from the regularity of the functions involved and their growth behavior. Since 
B + C + k I is bounded, the Lumcr-Phillips Theorem gives that the opera- 
tor A + B + C + kX generates a quasi-contraction semigroup. Specifically, the 
semigroup T = {T(t)} t >o, generated by the operator A + B + C, satisfies 

\\T(t)\\ < e~ Kt , t>0. (2.24) 

Finally, if condition (12. lip holds, we can choose K > 0. Hence the semigroup 
T = {T(t)}t>o is uniformly exponentially stable. 

□ 

Corollary 2.2 A stationary solution of Eqs. (|l.l[) - (ll.5|) is linearly asymp- 
totically stable if condition (|2.11|) holds true. 

Remark 2.3 For the stability of the trivial equilibrium n* = 0, the criterion 
(|2.11|) reduces to 

H(s,0) > P(s), se[0,oo). (2.25) 



Since 



m=i M<_['Kv,o)+o dy>d8 



< 



o t(s.o) I Jo 7(y,o) 

J /2(s,0) ex p|-^ /"M)^} ds = l (2.26) 
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with /2(s,0) = ^(s,0)/7(s,0), (j2T2"5j) implies i?(0) < 1. This is the well-known 
stability criterion of the trivial steady state in scramble competition. In the 
case of the hierarchical model discussed in [T3] our stability condition for the 
trivial steady state implied i?(0) < 1 as well. 



In this section we will characterize part of the point spectrum of the linear 
semigroup generator when the attack rate a assumes a special form. This result 
allows us to establish an instability result, thus complementing our stability 
result in the previous section. 

Throughout this section we make the assumption that the attack rate is 
separable, i.e. 



We can interpret a 2 (x 2 ) as a measure for the likelihood that individuals of size 
x 2 attack, while ai{x\) represents the likelihood of being attacked at size x\. 
Note that condition (|3.1|) assumes no correlation between these two events. Even 
though this assumption might appear as unsatisfactory from a biological point 
of view, it makes, however, analytical progress possible and henceforth has the 
potential to shed light on the case of a general attack rate. The main difficulty 
in case of general attack rates is that the operator C need not be a finite rank 
operator. However, our characterization of the point spectrum of the linearized 
operator A + B + C in the following developments relies essentially on this fact. 

We also make the biologically plausible assumption that there is a constant 
/io > such that 



Clearly, this assumption is automatically guaranteed if we suppose that \x > /io- 
The particular choice of the attack rate (|3.1[) allows us to cast the operator 
C in the form 



3. Instability via eigenvalues 



a(xi,x 2 ) = a 1 (xi) a 2 (x 2 ) 



(xi,x 2 ) G [0, oo ) x [0, oo). 



(3.1) 



(jl(s,E*(s)) +M*(s) > fi for all s £ [0, oo). 



(3.2) 



Cu = -uigi - u 2 g 2 , 



(3.3) 



where we define 




g 2 = ax n*. 



(3.4) 



Hence the eigenvalue problem 



(A + B + C- XX)u = 



u(0) = A(u) 
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assumes the form 

7* (s) u s + 0* (s) + A) it(s) + gi (s) Mi + g 2 (s) C 2 = 0, (3.5) 

= / "T77vT M ( s ) ds -53«i, (3-6) 
Jo 7 (0) 

with 

33 =7s(0,£*(0)K(0)a 2 (0). 

For 

ReA>-^ , (3.7) 
the solution of equation (|3.5p is found to be 

/* s (j (r) 

u(s) = u(0) 7r(s, A) - tti tt(s, A) / / 1 — -^dr 

Jo J*{r)n(r,X) 

-u 2 n(s,\) [ S g2(r) dr, (3.8) 
7 7*(r)7r(r,A) 

where we have made use of the notation 

r V0/) + A 



?r(s, A) = cxpj-^ 



7*(y) 



di/ 



7*(0) exp f pW)) + ^)+A (39) 



7*0) I Jo i*{y) 

Note that condition (|3.7p ensures that u, given by Eq. (|3.8[) . is in M /1,1 (0,oo). 
We multiply equation (|3.8[) by c{s)a^{s) and 0:2 (s)j respectively, and integrate 
from zero to infinity to obtain 

u(0) ai(A) + ui(l + 02(A)) + u 2 aa(A) - 0, (3.10) 
u(0) 04(A) + fli 05(A) + m 2 (1 + oe(A)) = 0, (3.11) 

where 

POO 

ai(A) = — / c(s) ai(s) 7r(s, A) ds, 
Jo 

fl2(A) = / C(S)Q1(S) " (S ' A) I r W^,A) drds ' 

f fla(r) 
7*W7r(r,A) 



03(A) = / c(s) ai(s) 7r(s, A) / TTT~T w drds > 



04(A) = — / 0:2(3) 7r(s, A) ds, 



05(A) = / a 2 (s)7r(s,A) / —^^-—drds, 

f ga(r) 

7*('") 7r (»",A) 



06(A) = / a 2 (s) 7r(s, A) / — ^ 2 - j — —drds. 
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Next we insert the solution (|3.8[) into the boundary condition (|3.6[) to obtain 
u(0) (1 + 07(A)) + Si MA) + .93) + u 2 09(A) = 0, (3.12) 

where 

07(A) =-y -^y7r(s,A)ds, 

a8(A) = i r¥) <sA) k r(^M Ws ' 

a 9 ( A )= rM^n^X) f^#-W, 

Now we can give a characterization of part of the point spectrum of the operator 
A + B + C. 

Theorem 3.1 For any A G {z 6 C | Rcz > — /io}j we Ziaue A £ <7 P (.A + 6 + C) 
if and only if A satisfies the equation 

( ai (A) 1 + 02(A) 03(A) \ 
X(A) = dot 04(A) 05(A) l + a 6 (A) =0. (3.13) 

\ 1 + 07(A) g 3 + a 8 (A) 09(A) / 

Proof. If X E o- p (A + B + C) , then Eqs. (|3.5|) . (|3.6|) admit a nontrivial solution 
u. Hence for this A there exists a nonzero solution vector (w(0), ui, U2) such 
that Eqs. 1(330 1) - ((332 ) hold true. Thus if (A) = 0. Conversely, if K(X) = 
for some A and (u(0), u\, U2) is a nonzero solution of Eqs. (|3.10[) - (|3.12p . then u. 
given by Eq. (|3.8p . is a nonzero solution of Eqs. (|3.5p - (|3.6|) at least if u(0) ^ 0. 
If, however, u(0) — 0, the only possible scenario for u to vanish when defined 
by Eq. (|3.8| would be the condition that 

u 1 g 1 = -u 2 g2 (3-14) 

hold true. Then Eqs. (|3.10|) . ()3.1ip would immediately give u\ = = u-i in 
contradiction to our assumption on (u(0), u\, S2). 

□ 

Since 

/0 1 0\ 

lim if(A) = det 1=1, (3.15) 

the limit being taken in R, we can formulate the following simple instability 
criterion, which follows immediately from the Intermediate Value Theorem. 

Theorem 3.2 A stationary solution n* of Eqs. (|l.l[) - (jl.5[) is linearly unstable 
ifK(0)<0. 



Note that for n* = 

#(0) = 1 -i?(0). 



(3.16) 
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Hence the stationary solution = is linearly unstable if R(0) > 1. In the 
remainder of this section we will concentrate on positive stationary solutions. 

Example 3.3 Let us assume that the rate of an individual of size s to attack 
another individual is proportional to the product of the probability for an in- 
dividual of size s to be attacked and its energetic value. Mathematically, this 
condition is modeled by the relation 

c(s)«i(s) =pa 2 (s), se[0,oo), p£l + . (3.17) 

The constant p denotes here the proportionality factor. This condition is bi- 
ologically relevant since environmental pressure conceivably makes individuals 
of higher energetic value, which are usually larger in size, not only more prone 
to be attacked, but also more aggressive. In addition we will assume that the 
attack rate of minimal size individuals (newborns) equals zero, i.e. 0:2(0) = 0. 
With 

tt(s) d = 7r(s, 0) (3.18) 

we obtain 

K(Q)= a 2 (s)n(s)ds / tt(s) / ————drds. (3.19 

Jo Jo 7(0) 7o 7T)T( r ) 

Note that (3 ^ by Eq. (|1.8p . Hence the instability criterion of Theorem 13.21 is 
satisfied if, for s > 0, 



p 1 a 1 (s)+a' 2 (s)-f E (s,E*(s))+a 2 (s) (J^Je(s, E* (s)) 



where we have used the relation 

> s *( s ) + /i(s,£*( s )) + M*( s ; 



n*(s) = -n*(s) 



7*(s) 

This instability condition automatically excludes the case a 2 = 0. 

Remark 3.4 In the age-structured case where 7 = 1 and where we may use 
a for age in lieu of size s we find for the preceding example 

M(a,t) = E(a,t)p- 1 ^-, (3.21) 

assuming that a 2 does not vanish. Hence the net reproduction function R(n) 
- as in the case of scramble competition - can be considered a function of the 
environment E, i.e. 

R(n) = R(E) = J™ (3(a) exp {- £ M (y, E(y, •)) + E(y, •) p^ da . 

(3.22) 
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When interpreting R as a nonlinear operator between sets of bounded continu- 
ous functions and ignoring all issues pertaining to regularity, we may formally 
deduce the Frechet derivative of R at a stationary state E* . The result of this 
formal calculation is 

Re(E*)=- [ /3(a)exp(-/ rtv,E.(y)) + E.(y)p- 1 ?^dy 

Jo I Jo My) 

x (j\ E {y^{y))+p- l ^-dy\ da. (3.23) 

Since in the scenario considered the instability condition (|3.20[) reduces to 

p _1 ai(a) +a 2 (a)u, E (a, E*(a)) < 0, oe[0,oo), (3.24) 
Eq. (j3~2"3"j) implies 

R E (E*) > 0. (3.25) 

Hence in this special case of model ingredients condition (|3.20|1 allows a formal, 
but intuitively clear biological interpretation, similar to the scramble competi- 
tion case in [12] : if the net reproduction rate R is increasing at a stationary 
environment £7*, then the equilibrium is unstable. 

4. Further spectral analysis 

Throughout this section we consider a positive stationary solution and 
assume that condition (13.21) holds true. 



Theorem 4.1 

o(A + B) n {z e C | Re z > -no} = P U {0}, (4.1) 

where the set P U {0} consists of simple, isolated eigenvalues A of A + B such 
that 

- fi < Re A < for Xe P. (4.2) 



A(exp<j-/ £bl±± dr}) ±\. (4.3) 



Proof. Suppose first that AeC with Re A > — no is such that 

*(r) + 
7*0") 

Then the equation 

(XI-(A + B))u = f (4.4) 
with / S i 1 (0, oo ) has the unique solution 

«(«) = «(0) exp(- f P * {r } ( + , X dr 
I Jo 7 W 

(0 + 

7*W J 7*(y) 



^(-/■^W+A^/fe)^ (4 . 5) 



I Jy 
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where 



«(0) = 








^ exp ■ 


u 


, T{r) J 


1 /(») 

7*(y) 




!-A( 


exp ■ 


("J 


o 7*(r) J 








(4.6) 



Note that our condition on A ensures that A can be applied to deduce (|4.6|) . 
Moreover, it readily follows that u belongs to W 1 ' (0, oo). Hence we conclude 
that A G p(A + B). Now suppose that A G C with Re A > — fiQ is such that 

A (-H^})- 

In light of Eq. (|1.8|) . A = solves (|4.7|) . In fact, A = is the only solution with 
Re A > 0. Let 

^H-jC^*})- 1 <"» 

for Re A > — fiQ. Then L is analytic for Re A > —fiQ and 

Jo 7*(0) Wo 7*M J io 7*(0 

Hence is a pole of the resolvent operator of A + B. Any other zero of L 
gives rise to a pole of the resolvent operator of A + B since L is analytic and 
nonconstant in the simply connected set Re A > — jiQ. Finally, the representation 
of the resolvent operator, given through the solution u in (|4.5[) , (|4.6p , shows that 
the spectral projection associated with each pole of the resolvent operator is a 
rank one operator. 

□ 

When the operator C is compact as in the special class of attack rates given 
in (|3.3p . we immediately obtain the following result. 

Corollary 4.2 Suppose that the operator C is compact. Then a (A + B + C) D 
{z G C | Re z > — fig} consists of isolated eigenvalues ofA + B + C of finite 
multiplicity. 



5. Asynchronous exponential growth 

The purpose of this section is to gain deeper insight into asymptotic prop- 
erties of solutions of the linearized system (|2.ip - (|2.2[) . In particular, we are 
interested in solutions of the linearized problem which grow exponentially in 
time such that the proportion of individuals within any size range compared to 
the total population approaches a limiting value as time tends to infinity, inde- 
pendently of the size distribution of the initial population. This phenomenon is 
called asynchronous exponential growth and is known to be present, e.g., in the 
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age-structured case. Mathematical definitions will be given below. The prop- 
erty of asynchronous exponential growth is important insofar as solutions can be 
regarded as asymptotically factorizable (with respect to time and size). Popula- 
tions of this kind are often called ergodic (see pi]). We refer to [4] [T0 | [TB I \T§ \ [23] 
for this and related notions. 

In the framework of linear semigroup theory a strongly continuous semigroup 
5* = {S(t)} t>Q on a Banach space y with generator As and growth bound 

s (A s ) d = sup {Re A | A 6 u (As)} (5.1) 

is said to exhibit balanced exponential growth (BEG for short) if there exists a 
bounded linear projection n on y such that 

lim \\e- s{As) 1 S(t) - n|| = 0. (5.2) 

t — >oo 

The semigroup S = {S(t)} t>0 is said to exhibit asynchronous exponential 
growth (AEG for short) if it exhibits BEG with a rank one projection n. For 
positive semigroups there exist well-known characterizations of BEG and AEG, 
see [HE]]. Our analytical approach will be guided toward these results. 

Theorem 5.1 Suppose that, for every y > and a.e. s > 0, 

c (y) a (y, s) ((je(s, E*(s)) n*(s)) s + (j, E (s, E*(s)) n*(s)) 

+ c(y) D 2 a(y, s)j E (s, E*(s))n*(s) + a(s, y) n* (s) < 0, (5.3) 
^(s)-7 E (0, J B*(Q))n*(Q)c(s)a(s,0) > 0. (5.4) 

Then the semigroup {T(i)} t >o, generated by the operator A + B + C , is positive. 

Proof. Condition (|5.3p implies that the operator C is positive, hence we can 
restrict ourselves to the operator A + B. Condition ()5.4|) ensures that the func- 
tional A is nonnegative. Consequently, the solution (|4.5|) . (|4.6j) of the resolvent 
equation (|4.4|) is well-defined and nonnegative if A > is sufficiently large. 

□ 

Corollary 5.2 Suppose that condition (|5.4p holds true. Then the semigroup, 
generated by A + B, is positive. 

Let us recall a useful characterization of irreducibility on L 1 (f2, m) (see [TO]): 
A strongly continuous, positive semigroup S — {S(t)} t>Q on the Banach lattice 
y = L 1 (fi,m) with generator ^5 is irreducible if, for / e y with / > 0, 
(XI — As) 1 /( s ) > f° r Ti-almost all s € fl and some A > s (As) sufficiently 
large. 

Theorem 5.3 Suppose that the positivity conditions (|5.3|) . (|5.4p hold true. Then 
the semigroup T = {T(t)} t >o, generated by A + B + C, is irreducible. 
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Proof. Since C is positive and since, for A > sufficiently large, the resolvent 
operator of A + B is positive as a consequence of Corollary 15. 2i we deduce from 

oc 

(XI- (A + B + C))' 1 = ({^-{A + B)y l C) (Xl-(A + B)y 1 



(5.5) 

that it suffices to prove the irreducibility of the semigroup generated by A + B. 
This result, however, follows immediately from the representation of solutions 
of the resolvent equation ()4.4|) . given by (14. 5|) . (|4.6| . 

. □ 

Before we formulate the main result of this section, let us review the notions 
of essential norm, growth bound, and essential growth bound, and some of 
their properties. Our discussion follows closely [10]. Suppose that As is the 
infinitesimal generator of the strongly continuous semigroup S = {S(t)} t>0 on 
a Banach space [V- Then the growth bound of the semigroup is defined by 

def IgJIgCOjj ,,,, 

^o(As) = lim — : — • (5-6) 
For a bounded linear operator T on y, the essential norm is given by 

||T|| oss d ^ f dist (T, £(?)), (5.7) 

where IC(y) denotes the set of compact linear operators on y. Of course, the 
essential norm is generally not a norm on the set of bounded linear operators on 
y. It is, however, a norm on the Calkin algebra of y, see [10] and the references 
therein. Finally, the essential growth bound of the semigroup S = (S(t)) t>0 on 
y with generator As is defined by 

i a ^ dcf r HISWIIcss , K o\ 

io css (As) = hm . (5.8) 

t— >oo t 

It is readily seen that, for /C e JC(y), 

^css (As) (As + K) . (5.9) 

The significance of the essential growth bound lies in the central fact that 

uj q (As) = max {uj css (A s ) , s (A s )} ■ (5.10) 

Theorem 5.4 Given a positive stationary solution n*, suppose that conditions 
(|3.2p . (|5.3p . and (|5.4p ZioW trwe and t/iai the operator C is compact. If 

a (A + B + C) n {z e C\Rez > 0} ^ 0, (5.11) 

£/ien i/ie linear semigroup T = {T(t)}t>o, generated by A + B + C, exhibits 
AEG. 
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Proof. First we note that A + B has nonempty spectrum and generates a posi- 
tive semigroup by Theorem l4.1l and Corollarv l5.2l Hence Derndinger's Theorem 
(see [SlUn]) proves that 

u o (A+B) = s(A+B) = 0, (5.12) 

where the last equality follows from Theorem 14.11 . Similarly, we obtain from 
Theorem 15. II and Derndinger's Theorem that 

u (A + B + C) = s(A + B + C) >0. (5.13) 

Here the last inequality is given by assumption (|5.1ip . Consequently, in light of 
Eqs. (5J§, (GEH), we have 

lu css (A + B + C)= uj css (A + B)< lu {A + B) = < lj q (A + B + C). (5.14) 

Hence by Theorem 15. II and Theorem 15.31 the semigroup T = {T(t)} t >o is posi- 
tive and irreducible with essential growth bound strictly smaller than its growth 
bound. The claim follows now immediately from Theorems 9.10 and 9.11 in [4]. 

□ 

6. Conclusion 

In this work we have studied the asymptotic behavior of solutions of a lin- 
earized size-structured cannibalism model, recently introduced in [7j . The vital 
rates in this model depend on a structuring variable (size), which takes values 
in an unbounded set, and on an infinite dimensional interaction variable (en- 
vironment), describing the environmental feedback on individuals. Population 
models of this type are notoriously difficult to analyze. The reason for this 
difficulty is that the essential spectrum of the linearized operator is typically 
not empty, and even the point spectrum cannot be characterized in general via 
zeros of a characteristic function. The latter obstacle has already been observed 
in 13J for a similar quasilinear hierarchically size-structured model. Therefore, 
analytical results, in particular with respect to the qualitative behavior of so- 
lutions, are rather rare in the literature, at least to our knowledge. We would 
like to point out that the emphasis in the present work was to demonstrate how 
analytical techniques can be developed and used to treat qualitative questions 
of physiologically structured population models, where the structuring variable 
is unbounded and competition is incorporated through infinite dimensional in- 
teraction variables. 

Here, using two different strategies, we have formulated linear stability and 
instability criteria for equilibrium solutions of the model. We derived a general 
instability criterion in the case of a separable attack rate and extended our dissi- 
pativity approach, employed in [13) for a model with finite size span, to the case 
of infinite size span. We also carried out a more refined spectral analysis of the 
linearized operator which allowed us to gain deeper insights into the asymptotic 
behavior of solutions of the linearized system. In particular, we investigated the 
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question whether solutions of the linearized problem exhibit asynchronous ex- 
ponential growth and gave sufficient conditions for an affirmative answer. There 
are only few results of this type for models with unbounded structuring variable, 
see [TH [201 EU [22] . In passing, we also note that our spectral analysis gives rise 
to an extension of our stability results in |12j for a model with finite size span 
to one with infinite size span. 

Acknowledgment 

JZF was supported by the EPSRC grant EP/F025599/1. TH gratefully ac- 
knowledges support through NSF-Grant DMS 0709197. We thank the reviewers 
for their thoughtful comments and suggestions. 



References 

A. S. Ackleh and K. Ito, Measure- valued solutions for a hierarchically 
size-structured population J. Differential Equations 217 (2005), 431-455. 

R. Borges, A. Calsina and S. Cuadrado, Equilibria of a cyclin structured 
cell population model, Prebulicacions (2008), Universitat Autonoma de 
Barcelona. 

A. Calsina and J. Saldaha, Asymptotic behavior of a model of hierarchically 
structured population dynamics, J. Math. Biol. 35 (1997), 967-987. 

Ph. Clement, H. J. A. M Heijmans, S. Angenent, C. J. van Duijn, and B. de 
Pagter, One-Parameter Semigroups, North-Holland, Amsterdam 1987. 

J. M. Cushing, A size-structured model for cannibalism, Theoret. Popula- 
tion Biol. 42 (1992), 347-361. 

R. Derndinger, Uber das Spektrum positiver Generatoren, Math. Z. 172 
(1980), 281-293. 

O. Diekmann and M. Gyllenberg, Abstract delay equations inspired by 
population dynamics, in "Functional Analysis and Evolution Equations" 
(Eds. H. Amann, W. Arendt, M. Hieber, F. Neubrander, S. Nicaise and J. 
von Below), Birkhauser, (2007), 187-200. 

O. Diekmann, Ph. Getto, and M. Gyllenberg, Stability and bifurcation 
analysis of Volterra functional equations in the light of suns and stars, 
SIAM J. Math. Anal. 39 (2007), 1023-1069. 

O. Diekmann, M. Gyllenberg, and J. A. J. Metz, Steady state analysis of 
structured population models Theoret. Population Biol. 63 (2003), 309-338. 

K.-J. Engel and R. Nagel, One-Parameter Semigroups for Linear Evolution 
Equations, Springer, New York 2000. 



Asymptotic analysis of a size-structured cannibalism model 



19 



[11] J. Z. Farkas and T. Hagen, Stability and positivity results for a generalized 
size-structured Daphnia model with inflow, Appl. Anal. 86 (2007), 1087- 
1103. 

[12] J. Z. Farkas and T. Hagen, Stability and regularity results for a size- 
structured population model, J. Math. Anal. Appl. 328 (2007), 119-136. 

[13] J. Z. Farkas and T. Hagen, Hierarchical size-structured populations: The 
linearized semigroup approach, submitted. 

[14] Ph. Getto, O. Dickmann, and A. M. de Roos, On the (dis)advantages of 
cannibalism J. Math. Biol. 51 (2005), 695-712. 

[15] M. Gyllenberg and G. F. Webb, Asynchronous exponential growth of semi- 
groups of nonlinear operators, J. Math. Anal. Appl. 167 (1992), 443-467. 

[16] T. Hagen, On the semigroup of linearized forced elongation, Appl. Math. 
Lett. 18 (2005), 667-672. 

[17] M. Iannelli, Mathematical Theory of Age-Structured Population Dynamics, 
Giardini Editori, Pisa 1994. 

[18] A. Pazy, Semigroups of Linear Operators and Applications to Partial Dif- 
ferential Equations, Springer, New York 1983. 

[19] S. Piazzera and L. Tonetto, Asynchronous exponential growth for an age 
dependent population equation with delayed birth process, J. Evol. Equ. 
5 (2005), 61-77. 

[20] J. Priifi, Stability analysis for equilibria in age-specific population dynam- 
ics, Nonlin. Anal. 7 (1983), 1291-1313. 

[21] H. R. Thicme, Positive perturbation of operator semigroups: growth 
bounds, essential compactness, and asynchronous exponential growth, Dis- 
crete Contin. Dyn. Syst. 4 (1998), 735-764. 

[22] H. R. Thieme, Balanced exponential growth of operator semigroups, J. 
Math. Anal. Appl. 223 (1998), 30-49. 

[23] G.Webb, Theory of Nonlinear Age-Dependent Population Dynamics, Mar- 
cel Dekker, New York 1985. 



